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ABSTRACT 

We have applied so-called 'optimised' perturbation theory to resolve the 
renormalisation-scale (//) ambiguity of exact O(a^) QCD calculations of event 
shape observables in e"'"e~ hadrons. We fitted the optimised predictions for 15 
observables to hadronic Z*^ decay data from the SLD experiment to determine 
as{M^). Comparing with results using the physical scale fi = Mz we found no 
reduction in the scatter among aa[M'^) values from the 15 observables, implying that 
the 0(0;^) predictions with optimised scales are numerically no closer to the exact 
all-orders results than those with the physical scale. 
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1. Introduction 

The theory of strong interactions, Quantum Chromodynamics (QCD) [1], contains 
in principle only one free parameter, the strong interaction scale A. Tests of QCD 
hence comprise comparisons of measurements of A in different processes and at 
different hard scales Q. In practise QCD calculations of observables are performed 
using finite-order perturbation theory, and calculations beyond leading order depend 
on the renormalisation scheme employed, implying a scheme-dependent A. Here we 
consider calculations performed in the modified minimal subtraction scheme {MS 
scheme) [2], and use the strong interaction scale I^jjg for five active quark fiavours. 
If one knows one may calculate the strong coupling ag{Q'^) from the solution 

of the QCD renormalisation group equation [3]. Because of the large data samples 
taken in e"'"e~ annihilation at the resonance, it has become conventional to use 
as a yardstick aa{M'^), where Mz is the mass of the boson; Mz ~ 91.2 GeV [4]. 
Tests of QCD can therefore be quantified in terms of the consistency of the values of 
q;s(M|) measured in different experiments. 

Measurements of aa{M'^) have been performed in e"'"e~ annihilation, 
hadron-hadron collisions, and in deep-inelastic lepton-hadron scattering experiments, 
covering a range of from roughly 1 to 10^ GeV^; for recent reviews see [5,6]. Within 
the errors there is a remarkable degree of consistency between these measurements; 
an average yields q;s(M|) = 0.117 ± 0.006 [6], implying that QCD has been tested 
to a precision of about 5%. This precision is however rather modest compared with 
the achievement of sub-1% level tests of the electroweak theory [7], and this is due 
primarily to the theoretical uncertainties that dominate most of the experimental 
measurements. These uncertainties are due to both the restriction of complete 
perturbative QCD calculations to low order, and non-perturbative effects that are 
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presently incalculable in QCD. The latter are often called 'hadronisation uncertainties' 
or 'higher twist effects'. 

For any observable, truncation of the perturbation series at finite order causes 
a residual dependence on the (scheme-dependent) renormalisation scale ji. In the 
case of e+e" annihilation, perturbative QCD calculations of hadronic event shape 
observables, such as jet rates [8] or thrust [9], have been performed exactly only 
up to second order in ag [10,11], and have been used extensively by experiments 
at the PETRA, PEP, TRISTAN, SLC, and LEP colliders for measurement of 
q;s(M|) [5,6]. The precision of these measurements is limited predominantly by the 
lack of knowledge of higher-order contributions, whose effect can, by definition, only 
be estimated in an ad hoc fashion. A consensus has arisen among experimentalists 
to estimate this effect from the //-dependence of the ^jfg, or equivalently aa{M'^), 
values derived from fits of the calculations to the data, see eg. [12], and to quote a 
corresponding renormalisation scale uncertainty. This procedure, well-motivated in 
that the //-dependence caused by the truncation of the perturbation series would be 
cancelled by addition of the higher-order terms, is, however, arbitrary, and is not 
equivalent to knowledge of the size of the a priori unknown terms. This arbitrariness 
is manifested in the fact that the experimental collaborations have chosen different 
ranges over which to vary /t; combined with their different choices of observables 
and averaging methods, this has led to the variation among quoted central values 
of as(M|) and scale uncertainties shown in Table 1. Despite this variation it is 
apparent that the scale uncertainty is much larger than both the experimental error 
and the hadronisation uncertainty, and represents the most serious limitation towards 
improved precision on measurements of as{M^) using hadronic event shapes in e"'"e~ 
annihilation. 
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The renormalisation scheme/scale ambiguity of QCD calculations has been 
discussed extensively in the literature [17]. The scheme ambiguity can be 
sidestepped by simply adopting one particular scheme, usually MS, as a reference; a 
next-to-leading order calculation of an observable in any scheme can then always be 
translated to the reference scheme [18]. The best resolution of the scale ambiguity 
would be to reduce its effect by calculating observables to higher order in perturbation 
theory. Though this is in principle possible, the large number of Feynman diagrams 
involved renders the task difficult and unattractive. In e"'"e~ annihilation only two 
inclusive observables, the hadronic cross-section ratio R [19] and the r hadronic 
decay ratio Rr [20], have been calculated exactly up to 0(0;^). For the hadronic 
event shape observables O(tt^) contributions have not yet been calculated completely, 
although progress has been made recently for some observables in the form of 
all-orders 'resummation' of leading and next-to-leading logarithmic contributions in 
the two-jet ('Sudakov') region [21]. These calculations have been 'matched' to the 
exact O(a^) results to yield improved predictions which, though not formally complete 
at any order beyond the second, have been found to yield a reduced dependence 
on the renormalisation scale, but at the cost of additional uncertainties relating to 
ambiguities at 0(0;^) in the matching procedure; see eg. [12]. In this paper we 
consider alternative approaches which attempt to 'optimise' O(a^) QCD predictions 
by choosing the renormalisation scale according to ad hoc criteria. 

2. Optimised Perturbation Theory 

The O(a^) prediction for an infra-red- and coUinear-safe observable X can be 
written: 

ax = -^{X,f,)=Ws{i,)A{X) + a;Hfi){A{X)27rbo\n{fi''/Q^) + B{X,Nf)) (1) 
a aX 
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where = aa/2'K, ba = (33 — 2Nf)/12'K, Nf is the number of active quark 
flavours, Q is the cm. energy, and A[X) and B{X,Nf) can be calculated [11]. 
Explicit dependence on the renormalisation scale jj, can be seen in the next-to-leading 
coefficient. We consider five possibilities for choosing this scale: 

(i) Physical scale: ii = Q, the hard scale of the interaction; in e"'"e~ annihilation at 
the resonance this corresponds to ji = Mz- This choice explicitly removes the 
logarithmic term in Eq. (1). 

(ii) Experimentally-optimised scale: yu can be derived from a simultaneous fit of 
^Ws ^° experimental data. This is entirely pragmatic and allows the data to 
'choose' the scale. The resulting Aj^-g and values are highly correlated [22]. 

(iii) PMS scale: Since the all-orders result would be independent of 
renormalisation scale, Stevenson suggests that /j, be chosen according to the 'Principle 
of Minimal Sensitivity' (PMS) [23], from the solution of: 

If - »• 

(iv) FAC scale: Grunberg suggests that /j, be chosen to give the 'fastest apparent 
convergence' (FAC) of the series [24], so that the second-order term in eq. (1) vanishes: 

A(X) 27r6o ln(AtV<3^) + B{X,Nf ) = 0. (3) 

At next-to-leading order this is equivalent to the 'effective charge' (EC) approach 
[24,25]. 

(v) BLM scale: Brodsky, Lepage and Mackenzie advocate [26] that n be chosen 
to remove the A^y^-dependence of the second-order term in eq. (1): 



11 = Qexj){3{B{X,Nf = 5)-B{X,Nf = 4)/2A{X))} (4) 
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As in Quantum Electrodynamics, this effectively incorporates quark and gluon 
vacuum polarisation contributions into the definition of the strong coupling. 

The experimentally-optimised, PMS, FAC and BLM approaches (ii)-(v) are usually 
collectively termed 'optimised' perturbation theory. In the PMS, FAC and BLM cases 
the optimised scale implicitly depends upon the value of the observable X. 

Early theoretical studies applied optimised perturbation theory to jet rates in 
e"'"e~ — > hadrons, and employed the PMS approach at Q ~ 34 GeV [27], and 
the PMS, FAC and BLM approaches at Q = Mz [28]. Experimentally-optimised 
scales were also determined for jet rates at Q = 29 GeV [29]. The focus of 
these early studies was largely on obtaining an improved description of the rate of 
4-jet production by 0{a1) QCD, which, for the physical scale Q, had been shown 
not to reproduce the PETRA data [8]. Until the SLC/LEP era, the influence of 
variation of the renormalisation scale on ag{M'^) measurements was usually not 
considered, corresponding, de facto, to choice of the physical scale. Early SLC/LEP 
q;s(M|) measurements based on jet rates and energy-energy correlations included 
experimentally-optimised scales [30] as well as the PMS, FAC and BLM methods 
[31,32]. However, theoretical controversy (see eg. [33]) motivated experimental groups 
to avoid specific scale-choice prescriptions, and subsequently to adopt the pragmatic 
approach of quoting an uncertainty on as{M^) by varying the renormalisation scale 
over a wide range; this approach has itself been criticised by theoreticians [34]. We 
are not aware of any comprehensive application of optimised perturbation theory to 
hadronic event shape observables, and of its effect on ^^(MI) determinations. 

Here we present such a study. We have calculated the PMS, FAC and BLM scales 
for the 15 collinear- and infra-red-safe hadronic event shape observables considered 
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in the SLD q;s(M|) measurement [12]. Using these scales, as well as the physical and 
experimentally-optimised scales, we have extracted as{M^) values from comparison 
with the SLD data. For each scale choice we have studied the scatter among the 
as{M^) values from the 15 observables, which one expects a priori to be reduced if 
the optimised perturbation series up to O(a^) are indeed better approximations to 
the all-orders results. 

3. Hadronic Event Shape Observables 

First we review briefly the hadronic event shape observables. Thrust T is defined 

[9] 

r = max ^il^;?^' , (5) 

where pi is the momentum vector of particle i, and ny is the thrust axis to be 
determined. We define r = 1 — T. An be found to maximize the 

momentum sum transverse to ny, and defined to be perpendicular to 

the two axes ny and fimaj- The variables thrust-major T^aj and thrust-minor Tmin 
are obtained by replacing Ht in Eq. (5) by n^aj or fimin, respectively. The oblateness 
O is then given by [35] O = Tmaj—Tmin- The C-parameter, C = 3(AiA2+A2A3+A3Ai), 
is derived from the eigenvalues Aj (i = 1,2,3) of the infrared-safe momentum tensor 
[36]: 

~ V I r?- I ' 

where p^ is the p-th component of the three momentum of particle i, and i runs over 
all the final state particles. 

Events can be divided into two hemispheres, a and b, of invariant mass Ma and 
Mfj, by a plane perpendicular to the thrust axis ny. The heavy jet mass Mjj is then 



defined [37] Mjj = max(Ma, Mj). Here we consider the normalized quantity p = -^r^, 

vis 

where E^is is the visible energy measured in each hadronic event. Jet broadening 
measures have been proposed in Ref. [38]. In each hemisphere a, b: 

is calculated. The total jet broadening Bt and wide jet broadening Bw are defined 
by Bt = Ba + Bfj and Byy = max(i?a, Bj,), respectively. 

For back-to-back two-parton final states r, O, C, Bt and Bw are zero; for planar 
three-parton final states 0<t<1/3, 0<O< l/\/3 and < C < 2/3; spherical 
events have r = 1/2 and C = 1. 

Another useful method of classifying the structure of hadronic final states is in 
terms of jets. Jets may be reconstructed using iterative clustering algorithms in which 
a measure yij, such as scaled invariant mass, is calculated for all pairs of particles i and 
j, and the pair with the smallest yij is combined into a single particle. This procedure 
is repeated until all pairs have yij exceeding a value ycut, and the jet multiplicity of 
the event is defined as the number of particles remaining. The n-jet rate Rniycut) is 
the fraction of events classified as n-jet, and the differential 2-jet rate is defined [31] 

n / X _ R^iVcut) - R2{ycut - ^Vcut) 

^2[ycut) = T- ■ («j 

^Ucut 

Several schemes have been proposed comprising different yij definitions and 
recombination procedures. We have applied the E, EO, P, and PO variations of the 
JADE algorithm [8] as well as the Durham (D) and Geneva (G) schemes [39]. 

Hadronic events can also be classified in terms of inclusive two-particle 
correlations. The energy-energy correlation (EEC) [40] is the normalized 



energy-weighted cross section defined in terms of the angle Xij between two particles 
i and j in an event: 

events ^ tj "^'S 
X~ 2 

where x (0 < x ^ 180°) is an opening angle to be studied for the correlations, 
Ax is the angular bin width, and Ei and Ej are the energies of particles i and j 
respectively. The shape of the EEC in the central region, x ^ 90°, is determined 
by hard gluon emission. The asymmetry of the EEC (AEEC) is defined as 
AEEC{x) = EEC {180° - x) - EEC{x). Another procedure, related to the angle of 
particle emission, is to integrate the energy within a conical shell of opening angle x 
about the thrust axis. The Jet Cone Energy Fraction (JCEF) is defined [41]: 

JCef{x)^j^^J: / E#^^(x'-x.)dx', (10) 

l^events^X ^^^^^^ ■ ^vis 

X 2 

where Xi = arccos (pj ■ ny/ \ pi |) is the opening angle between a particle and the 
thrust axis vector, Ht, whose direction is defined to point from the heavy jet mass 
hemisphere to the light jet mass hemisphere. Hard gluon emission contributes to the 
region corresponding to the heavy jet mass hemisphere, 90° < x ^ 180°. 

4. Measurement of as(M|) 

Distributions of these 15 event shape observables were measured [12] using a 
sample of approximately 50,000 hadronic decay events collected by the SLD 
experiment. The data were corrected [12] for detector bias effects such as acceptance, 
resolution, and inefficiency, as well as for the effects of initial-state radiation and 
hadronisation, to arrive at 'parton-level' distributions. 
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For each observable we employed the EVENT program [42] to calculate the 
coefficients A and B in Eq. (1). We then fitted the 0{al) calculation to the measured 
parton-level distributions, using first the physical scale ji = Mz, by minimising 
w.r.t. variation of A-^^. In each case the fit range was chosen so as to exclude the 
2-jet region, where resummation [21] of higher-order perturbative contributions is 
required [12], as well as the 4-jet region, where the 0{a,^) calculation is not expected 
to reproduce the data accurately; these ranges are indicated in Fig. 1. Each resulting 
value was translated into Q,g(M^); these are shown, with experimental errors 
[12], in Fig. 2(a). It can be seen that there is considerable scatter among the 15 
q;s(M|) values. Since the same data sample was used to measure each observable, 
and since the observables are highly correlated, this scatter is very significant. Similar 
results have been observed previously [16]. The scatter can be interpreted as arising 
from uncalculated higher-order perturbative QCD contributions, which a priori may 
be of different sign and magnitude for the different observables. Taking an unweighted 
average over all 15 q;s(M|) values, and a corresponding r.m.s. deviation, yields^: 

q;s(M|) = 0.1265 ± 0.0076 (physical scale) . 

We repeated this procedure using the experimentally-optimised-scale approach. 
In this case a simultaneous fit of and /i to each distribution was performed; 

the fitted value of fi is indicated in Fig. 1. The resulting pairs of A-^^ and fi values 
were translated to ^^(MI), which are also shown in Fig. 2(a) . Again, there is large 
scatter among the 15 a,s(Af|) values. It should be noted that for most observables 

^ A weighted average based on experimental errors yields as(Mf ) = 0.1273, which agrees with 
the unweighted average by less than the statistical error on as{M^) from a single observable. 

^ For the observable no minimum in w.r.t. variation of /« in the range lO^** < /M^ < 
10^ could be found. 
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the experimentally-optimised scale yields a lower value of as(M|) than the physical 
scale; this is because the optimised scale is typically smaller than Mz, which usually 
requires a smaller value of in order to fit the data [22]. For each observable 
such a systematic difference is encompassed by the renormalisation scale uncertainty 
on q;s(M|) considered in [12], which is also shown in Fig. 2. Taking an unweighted 
average and r.m.s. deviation yields^: 

aa{Mz) = 0.1173 ±0.0071 (experimentally — optimised scale). 

As expected from the preceeding discussion, the central value is lower than 
for the physical scale choice. However, the r.m.s. deviation obtained with 
experimentally-optimised scales is comparable with that resulting from choice of the 
physical scale. 

For each observable we then calculated the PMS, FAC and BLM optimised scales 
by solving Eqs. (2), (3) and (4), respectively; in the BLM case the next-to-leading 
coefficients B were calculated separately using the EVENT program for Nf = 4 and 
5. These scales are shown in Fig. 1. For each observable the following points are 
apparent: 1) the optimised scale depends strongly on the value of the observable; 2) 
with the exception of O, across the observable range the optimised scales are typically 
much smaller than the physical scale fi = Mz] 3) the PMS and FAC scales are almost 
identical; 4) as one approaches the 2-jet region, corresponding to X/{Xmax — Xmin) 
^ in Fig. 1(a), ycut ^ in Fig. 1(b), and x ^ 180° in Fig. 1(c), the BLM 
scale decreases whereas the PMS and FAC scales increase. In the last case the 
BLM behaviour conforms to the naive expectation that the optimised scale should be 

^ A weighted average using experimental errors yields as(M|) = 0.1166. 
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closely related to the momentum-transfer involved in the physical process, namely the 
radiation of soft and/or collinear gluons; the behaviour of the PMS and FAC scales 
in the two-jet limit does not appear to satisfy this expectation. 

For each observable we then fitted the O(a^) calculation to the measured 
distribution, using in turn the PMS, FAC and BLM scales, to determine A-^g^ and 
hence as(M|). The results are shown in Fig. 2(b); in the case of oblateness an 
acceptable fit with the BLM scale could not be obtained. For each observable 
it can be seen that the PMS- and FAC-derived Q!s(M|) values are very similar, 
whereas, typically, the BLM-derived as(M|) value differs from them. This behaviour 
follows from the correlation between the scale value (Fig. 1) and the corresponding 
^Ws I'squired to fit the data [22]. Comparing Figs. 2(a) and 2(b) it can be seen 
that for a given observable the PMS- and FAC-derived q;s(M|) values are often, 
though not always, close to that determined using the experimentally-optimised scale. 
Furthermore, for most observables the PMS-, FAC- and BLM-derived q;s(M|) values 
all lie within the range encompassed by the /i- variation considered in [12], though for 
p, Bw, D^, D[^, and (Bt), the BLM- (PMS/FAC-) derived values lie below this 
range. 

The most striking feature of Fig. 2(b) is that, for any of the PMS, FAC or BLM 
scale choices, there is considerable scatter among the ^^(MI) values from all the 
observables. In each case, taking an unweighted average over all the as(M|) values 
and a corresponding r.m.s. deviation yields^: 

a,(M|) = 0.1123 ±0.0079 (PMS scale) 

^ Weighted averages based on experimental errors yield central Q!s(Mf ) values of 0.1120 (PMS), 
0.1120 (FAC) and 0.1086 (BLM). 
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a,(M|) = 0.1123 ±0.0080 (FAC scale) 



a,(M|) = 0.1088 ±0.0075 (BLM scale). 

In each case the r.m.s. deviation is comparable with that resulting from choice of the 
physical scale, or of the experimentally-optimised scale. 

5. Summary and Conclusions 

We have determined q;s(M|) by fitting 0{a^) QCD predictions of 15 hadronic 
event shape observables to e"'"e~ annihilation data at the Z*^ resonance collected by 
the SLD experiment. We used five prescriptions for resolving the renormalisation scale 
ambiguity that arises in the truncated perturbative calculation, namely the physical, 
experimentally-optimised, PMS-, FAC- and BLM-optimised scales. The average 
q;s(M|) value, taken over all the observables, differs among these five procedures, 
which can be understood from the correlation of q;s(M|) with the renormalisation 
scale value [22]. More importantly, the scatter among the q;s(M|) values is equally 
large for all five prescriptions, the r.m.s. deviation being about 0.008. 

We conclude that the optimised perturbation theory approach does not reduce 
the scatter among the Cfes(M|) values determined from different observables. If such 
scatter is interpreted as arising from the effects of the uncalculated higher-order 
perturbative QCD contributions, then in the cases we have investigated these 
contributions appear to be as large for optimised scales as for the physical 
scale. Therefore, notwithstanding the possible merits of optimised perturbation 
theory from a theoretical point-of-view, we have demonstrated that this approach 
does not appear to offer any numerical advantage in terms of the accuracy of 
perturbative QCD predictions of e+e" event shapes. This is in agreement with the 
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expectations of a recent study by Barclay and Maxwell [18], who advocate the use 
of renormalisation-scheme-invariant quantities as probes of the size of uncalculated 
higher-order QCD effects. 

We thank our colleagues in the SLD Collaboration for support for this analysis. We 
also thank S. Brodsky, L. Dixon and C. Maxwell for helpful discussions and for their 
encouragement of this work. 
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Experiment 


Observables 


«,(M|) 


Exp. 


Errors 
Had. 


Scale 


Reference 


SLD 


T, 0, C, Ml Bt, Bw, 
D^, EEC, AEEC, JCEF 


0.123 


±0.003 


±0.002 


±0.011 


[12] 


ALEPH 




0.121 


±0.004 


±0.007 a) 


+0.007 
-0.012 


[13] 


DELPHI 


T, 0, C, Ml Ml , 
EEC, AEEC 


0.113 


±0.002 


±0.003 


±0.006 


[14] 


L3 b) 


i?f , EEC, AEEC 


0.118 


±0.004 


±0.004 


±0.006 


[15] 


OPAL 


r, 0, C, Mh, Md, m\^\ 

^^^d ^ ^2 ^ ^2 ' -'^2 ' 

L>f , AEEC, PTEC 


0.122 


±0.002 c) 




+0.006 ,^\ 
-0.005 ^) 


[16] 



Table l.as(M|) and errors from O(a^) QCD fits to hadronic event shape observables 
in decays. For each experiment results are taken from the most recent 

Qfg determination based on event shapes using O(a^) calculations. Definitions of the 
observables listed in the second column can be found in the references shown in the last 
column. The fourth column shows experimental errors, the fifth column hadronisation 
uncertainties and the sixth column scale uncertainties. Notes: a) uncertainty due 
to combined 'higher orders and hadronisation effects' [13], not including the scale 
uncertainty; b) we averaged the separate L3 measurements from jet rates and the EEC 
and AEEC; c) we estimated this value from Table 4 of Ref. [16]; d) Ref. [16] quotes a 
total uncertainty of ^Q QQg based on a weighted average over all 13 observables, taking 
correlations into account; subtracting the estimated experimental error in quadrature 
yields a theoretical uncertainty of ^Q QQg, which includes both hadronisation and scale 
uncertainties. 
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Figure Captions 



FIG. 1. Optimised values of the renormalisation scale plotted as / = ji'^/s. 
Experimentally-optimised scale (horizontal dotted line); PMS and FAC (solid), and 
BLM (dashed) scales. Differences between the PMS and FAC scales cannot be 
resolved in this figure. In (a) the scale is plotted vs. the dimensionless variable 
X/{Xmax - Xmin), where a: = r, p, Bt, Bw, O or C, and Xmax (Xmin) is the 
maximum (minimum) kinematically-allowed value of X. In (b) the scale is plotted 
vs. Hcut for D2{yc) calculated with the E, EO, P, PO, D and G jet-finding schemes. 
In (c) the scale is plotted vs. x for the EEC, AEEC and JCEF (see text). The 
range of each observable used in the fits to determine q;s(M|) is indicated by vertical 
dashed lines. 

FIG. 2. Values of tt,(M|) from QCD fits to the data using: (a) physical (solid 
circles), and experimentally-optimised (open circles) scales; (b) PMS- (solid circles), 
FAC- (solid triangles), and BLM- (open squares) optimised scales. In all cases only 
experimental error bars are shown. For each observable the shaded region indicates 
the total uncertainty estimated in Ref. [12], dominated by the contribution from wide 
variation of the renormalisation scale. 
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